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Abstract — We estimate the temperature distribution in the 
rails of an electromagnetic rail gun (EMG) due to the confinement 
of the current in a narrow surface layer resulting from the skin 
effect. In order to obtain analytic results, we assume a simple 
geometry for the rails, an electromagnetic skin effect boundary 
edge that propagates with the accelerating armature, and a 
current carrying channel controlled by magnetic field diffusion 
into the rails. We compute the temperature distribution in the 
rails at the time that the armature leaves the rails. For the range 
of exit velocities, from 1500 m/s to 5000 m/s, we find the highest 
temperatures are near the gun breech. After a single gun firing, 
the temperature reaches the melting temperature of the metal 
rails in a layer of finite thickness near the surface of the rails, 
for rails made of copper or tantalum. We plot the thickness of the 
melt layer as a function of position along the rails. In all cases, the 
thickness of the melt layer increases with gun velocity, making 
damage to the gun rails more likely at higher velocity. We also 
calculate the efficiency of the EMG as a function of gun velocity 
and find that the efficiency increases with increasing velocity, if 
the length of the gun is sufficiently long. The thickness of the 
melted layer also decreases with increasing rail length. Therefore, 
there is a tradeoff: for rails of sufficient length, the gun efficiency 
increases with increasing velocity but the melted layer thickness 
in the rails also increases. 

I. Introduction 

ELECTROMAGNETIC launch systems, such as the rail- 
gun, are based on transient phenomena [1]. During 
launch, the transient involves the build-up and penetration 
of a magnetic field into the surrounding metallic material. 
The dynamics of magnetic field penetration into the metal 
rails is described by a well-known diffusion equation [2]. The 
diffusion of the magnetic field leads to a skin effect, where 
a large current is transported inside a narrow channel. In a 
railgun that has a moving conducting armature, the effect is 
called a velocity skin effect (VSE) and is believed to be one of 
the major problems in limiting railgun performance [2], [3], 
because it leads to intense Joule heating of the conducting 
materials, such as rails and armatures. To what extent the 
VSE effect is responsible for limiting the performance of solid 
armatures is still the subject of research [4], [5]. 

Recently, using a new generation of magnetic field sensors, 
the magnetic field distributions caused by the VSE in the rails 
have been measured [6], [7]. These experimental efforts are 
even more significant in light of the large investments planned 
by the navy to develop EMGs and power sources for nuclear 
and conventional warships [8], [9]. The motivation for this 
paper is the large investment planned for EMG technology 



and the historical lack of understanding of the reasons for the 
low endurance of the gun rails in service, sometimes limited to 
one shot at maximum energies before replacement is needed. 

For high-performance EMGs, in order to increase the arma- 
ture velocity while keeping the length of the rails fixed, the 
current pulse during firing must be shorter and have a higher 
average amplitude, causing a stronger skin effect in the rails, 
which leads to an increase in Joule heating of the rails. In this 
paper, we show that the EMG efficiency is higher at higher 
velocity, but there is increased melting of the rails, leading to 
a tradeoff between efficiency of the EMG and melting of the 
rails due to Joule heating (during a single firing). We do not 
discuss gun barrel erosion due to repeated firings [10], [11]. 
We also do not treat the interaction of the high temperature 
plasma in the contact regions of the EMRG 

The article is organized as follows. We define our model for 
EMG heating in Section II. We consider the problem of rail 
heating in two steps. First, in Section III we assume a current- 
carrying channel described by a local skin depth S (along the 
rail) that has a simple time dependence due to motion of the 
armature, which leads to a time and position-dependent current 
density that causes the Joule heating. In Section IV, we use 
an improved expression for the current density based on the 
diffusion of the magnetic field into the rails. In Section V, we 
give a crude approximation to the armature heating. In Section 
VI, we discuss the efficiency of the EMG and its dependence 
on gun velocity and length. We show that the there is a tradeoff 
between efficiency and melting the rails. Finally, we present 
our conclusions in Section VII. 

II. Temperature Distribution in Rails 

The dynamics of an EMG can be properly described in 
terms of a thermodynamic free energy that expresses the 
coupling of the mechanical and electromagnetic degrees of 
freedom [12]. For example, for the case of an electromagnetic 
gun that can be described in terms of a lumped circuit model, 
which has a rotor coil with self inductance L R carrying current 
I R , and a stator coil, rail, and armature circuit with self 
inductance L with current /, the free energy has the form [13] 

F(I R , I, x, 9) = F (T) + \l r I 2 r + \l{x)I 2 + M{6) I R I 

(1) 

where L(x) depends on the a; -position of the armature, 

L(x) = L + L' x (2) 
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and where L' = dLgjdx, and Lq is the self inductance of the 
stator, rail and armature circuit when the armature is at x = 0. 
The interaction of the stator and rotor circuits is specified in 
terms of their mutual inductance, M(9), where 9 is the angle 
of the rotor coil with respect to the stator coil. The term F b (T) 
depends only on temperature. Derivatives of the free energy 
with respect to the coordinates, x and 0, give the generalized 
forces on the system [12]. By Newton's law, the acceleration 
of the armature, x, is given in terms of the derivative of the 
free energy with respect to the coordinate, 



dx J 



(3) 



where m is the mass of the armature (including the payload, 
comprising the launch package), leading to the well-known 
dynamical equation 

1 



m x 



I T 2 



L'l 



(4) 



Equation Q shows that the acceleration is directly propor- 
tional to the square of the instantaneous current, /. As a 
first approximation, we assume that the current in the rails 
is constant during a shot, I(t) = I a . If the current is constant, 
and the armature starts at t = at x = 0, and moves to the 
end of the gun rails at x = £ at time t — tf, then we have 
the following relations between armature mass m, gun length 
£, self inductance per unit length of armature travel L', and 
mass velocity v, 



Io = 

tf = 
t Q (x) 



_ _2_ ^zroj 1 / 2 



(5) 



The function t (x) gives the time at which the mass m is at 
position x along the rails, see section III. 

In order to calculate the Joule heating in the rails, a 
lumped circuit model is not sufficiently detailed. Instead, we 
must use a more detailed model, where the free energy is 
expressed in terms of electromagnetic fields with a spatial 
distribution. The problem is complicated due to the coupling of 
the electromagnetic and mechanical degrees of freedom. The 
problem is further complicated by the fact that the dynamics of 
an EMG shot is a transient effect in time. Consequently, when 
Long [14] and Nearing and Huerta [15] computed current 
density (and the heating in EMG rails) they assumed that 
the mechanical and electromagnetic degrees of freedom are 
decoupled. Furthermore, they assumed the armature was mov- 
ing at a constant speed so the problem became translationally 
invariant in time, thereby avoiding the complexities associated 
with the initial conditions and the resulting transient effects. 
In particular, it is the transient nature of the EMG shot that 
gives rise to a dependence of EMG performance on rail length, 



see for example our Eq. (14i and (36i for dependence of 



rail temperature rise and EMG efficiency, respectively, on 
rail length £. Using the above stated assumptions, Long [14] 
and Nearing and Huerta [15] solved for the complicated 
distribution of the current density using a simplified geometric 
model of the rails and armature [16]. 



In this paper, we use a simpler approach that addresses 
the transient nature of the EMG shot and allows us to get 
approximate analytic results for the temperature distribution in 
the rails and armature of an EMG. We model the skin effect 
in the rails, which limits the channel through which current 
can flow. For a given total current in the rails, a narrower 
channel (smaller skin depth) leads to a higher current density 
and results in greater Joule heating of the rails. In contrast to 
Nearing and Huerta [15], our approach allows us to discuss 
transient effects dependent on the length of EMG rails and how 



they impact EMG performance, see for example our Eq. ( 14 1 
and (36) . We use the same 2-dimensional simplified geometry 
for the rails and armature as Nearing and Huerta. However, 
we assume that the conducting armature is arbitrarily thin, 
which allows us to get simple results. Corrections to such an 
assumption are expected to be of order 0(a/£), where a is 
the armature width and £ is the length of the rails. We take 
the coordinate x running down the length of the rails, and we 
assume that the top rail occupies y > and the bottom rail 
occupies y < —6, where b is the rail separation. We take the 
length of the gun from breech to muzzle to be £. We assume 
that the rails have an arbitrary large thickness w in the z- 
direction, see Figure 1. When an EMG is fired, the current 
flows down one rail, through the conducting armature, and up 
the other rail to complete the electrical circuit. The transient 
response leads to a complicated distribution of eddy currents 
in the rails and the armature. We assume that we can describe 
this effect by the local skin effect that depends on position 
and time. 

We take the temperature rise due to the Joule heating to be 



T(x,y,t f ) - T = 



p C 



J 2 (x,y,t) 



dt 



(6) 



where T(x,y,tf) is the temperature at position (x,y) at time 
tf, when the armature leaves the rails, assuming the current 
starts at t = 0. The quantity To = T(x, y, 0) is the initial 
temperature at position (x,y) at time t = before the shot, 
C is the specific heat (assumed constant up to the melting 
point) of the metal rails, and p is the density of the metal 
rail. We assume the electric field E(x,y,t) is linearly related 
to the current density, J(x,y,t) = a E(x,y,t), where a is 
the electrical conductivity that is independent of temperature 
up to the melting point. In Eq. |6]), we have neglected the 
heat of melting, so the temperature rise is only valid up to 
the melting point of the metal rail. If latent heat of melting 
Lq is included, then the term —Lq/C must be added to the 
right side of Eq. ([6]). This additional term subtracts from the 
temperature rise that may be expected when additional Joule 
heat is created beyond what is required to reach the melting 
temperature. For Cu or Ta, this term is significant, with value 
Lq/C = 465 °C and 1141 °C, respectively. Below, we do not 
consider the temperature rise above the melting point of the 
metal rails. 
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III. Current Channel with Constant Current in 
Rails 

The moving armature carries all the current of the rails. As 
the armature moves, it exposes a new plane region on the rail 
that carries current. Due to the diffusive nature of the magnetic 
field H, and the relation 



J = curl H 



(7) 



current does not flow uniformly in the rail [12], [17], instead 
the current flows through a layer of thickness 5 that increases 
with time i, starting with 5 = as the leading edge of the 
armature passes a point on the rail. Assuming a plane geometry 
for the rail, we can approximate the time dependence of the 
skin depth to be [17] 



5(t) = 



At 



1/2 



(8) 



where /i is the magnetic permeability and a is the electrical 
conductivity, and t is the time elapsed since the armature has 
passed a given element on the rail. The exact factor inside the 
square root (here we take it to be 4) is somewhat arbitrary in 
defining a skin depth. (In Section IV we remove this arbitrari- 
ness by using the solution of the time-dependent magnetic field 
diffusion equation to compute the current density distribution.) 
As stated in the second to last paragraph in Section II, we 
assume an idealized armature that is arbitrarily thin in the 
^-direction and we assume that the skin depth 5 in the rail 
starts at zero thickness at the position where the armature 
contacts the rails. As remarked above, this drastic assumption 
is expected to have corrections of order 0(a/£) where a is 
the width of the armature and t is the length of the rail. This 
assumption allows us to obtain analytic results and see the 
dependence on a number of parameters. As a start, we model 
the conducting channel by assuming that the current density in 
the rail inside the skin depth S is a function of x and t but not 
y, and that the current density is zero outside the skin depth S, 
see Figure 1 [18]. Furthermore, we assume that the skin depth 
has zero thickness on the leading edge of the armature, since 
all the current in the rails has to flow through the armature. 
Assuming a current carrying channel of finite width 5, at time t 
and position x along the rail, the y coordinate of the boundary 
of the current carrying channel is 



Vc(x,t) 



S(t-t (x)) 
0, 



t > t (x) 
t < t (x) 



(9) 



where t (x) is the time the leading edge of the armature passes 
the position x. In other words, 



dt (x) 
dx 



1 



v{x) 



(10) 



where v(x) is the velocity of the armature when it is at point 
x, see Figure ([TJ. 

Consider an element of volume dV — wdxdy at position 
(x,y). Define r(x,y) as the time at which the boundary of 
the current carrying channel intersects this volume element. 



I 



x=0 
t= 



Fig. 1 . The rails and conducting armature are shown for the electromagnetic 
gun. At time t = the leading edge of the armature is at x = 0. The trailing 
edge of the armature leaves the rails at time t = tt. 




//////// 



< t < t{x, y) 



Fig. 2. At time t < r the skin depth boundary has not yet reached the 
element of volume dV = wdxdy at position (x, y). 



From Eq. Q for the time-dependent boundary of the current 
carrying channel, we find 



r(x,y) = t (x) + -^xay 



(11) 



For time < t < r(x,y) there is no current flowing through 
this volume element. At time t = r(x,y), the boundary of 
the current carrying channel intersects the volume element at 
position (x, y) and Joule heating starts. See Figures 2, 3 and 
4. Finally, at time t = tf the armature leaves the rails, the 
circuit is broken, and there is no more Joule heating of the 
element dV. 

The current density in the rail can be written down by 
considering three domain regions, see Figure 5. At position 
(x, y) and time t, we take the current density to have the form 

J(x,y,t) = 

x < x (t) and y < S(t — t {x)) 
x < x (t) and y > S(t — t a (x)) 
x > x (t) 

(12) 

where x a (t) is the function that gives the x coordinate of the 
armature at time t, and S(t) is given by Eq.([8ji. 

Using Eq. (|6|l and assuming a constant current in the 
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J = ! 




X 



x = x o (t) 



Fig. 3. At time t = r the skin depth boundary overlaps the element of 
volume dV = wdxdy at position (x, y). At this time, current starts to flow 
in the element dV and the temperature starts to rise due to Joule heating. 




t>T(x,y) 



Fig. 4. After the skin depth boundary passes the element of volume dV, 
heating of the element continues until the armature exits the rails at the time 
tf. For time t > tf, we assume that no energy is input into the element dV. 



armature, I(t) = Iq, from Eq. d6j) and (12i we obtain 



T(x,y,tf) - T Q = 



Aw 2 pC 



Log 



-'.(') 



-(x,y)-t a (x) 

0. 



y < 
y > 



— (*/ 

LLO \ J 
fl(T v J 



t (x)) 
t (x)) 



1/2 
1/2 

a 3) 

Equation (13i gives the temperature rise, T(x,y,tf) — To, 
of an element of volume at position (x, y) at time tf at which 
the armature leaves the rails, for a constant gun current over 
the time interval < t < tf. For realistic time-dependent 
currents, see Figure 3 in McCorkle [19]. 

Using a lump circuit model to describe the constant ac- 
celeration x = a of the armature for a constant current Iq, 
we can write the position of the armature as a function of 
time as x = hat 2 . The time at which the armature leaves the 
rails, £/, is then related to the length of the gun, £ = \at 2 . 
The time t a (x) at which the armature passes coordinate x 
is then t (x) = (2xja) Y l 2 , see also Eq. B. Using these 
approximations in Eq. ( [T3] l gives the temperature rise at 
position (x, y) 



T (x,y,tf) - T Q 



Log 



8 t 



0. 



4w 2 pC 



,1/2 



y < y(x) 
y > y(x) 



(14) 



Fig. 5. The three domain regions are shown for the current density. The non- 
zero current density is assumed to be inside the skin depth for x < x (t) 
where x (t) gives the x-coordinate of the armature at time t. 
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Fig. 6. Plot of Eq. [Bj, giving the temperature rise, assuming rails made 
from copper, plotted vs. position x along the rail. Parameters used are given 
in Tables I and II. The gun velocity is taken as 3000 m/s and the inital 
temperature T = 0°C. Constant current is assumed. 



where 



y(x) = 



1/2 



1- ~ 



(?) 



1/2 



1/2 



(15) 
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and x is assumed to be in the interval < x < £. Note that 
Eq. ( p"4} predicts that the temperature on the surface of the 
rail, at y = 0, is infinite. This feature of the solution is well- 
known and is not a problem [17]. The singularity with respect 
to y is integrable, and therefore the energy deposited in a thin 
layer near the surface is finite. We use Eq. ( 14 1 to locate the 



points (x, y) of the surface that reaches the melting point. 

Using the values for the gun parameters and material 
parameters in Tables I and II, we plot the distribution of the 
temperature rise in the rail, T(x, y,tf) — T , as a function of 
the y— coordinate, for several points along the rail given by 
coordinate x. For comparison, we also plot the melting point 
of Cu and Ta, which is 1083°C and 2996°C, respectively. 
The y— coordinate where the curve intersects the melting point 
of each metal is the depth to which melting of the rails 
occurs, if the EMG gun shot occurred with rails at initial 
temperature To- As mentioned previously, we have not taken 
into account the heat of melting, so these curves are not correct 
for temperatures above the melting point of the rails. 

Figure 6 and 7 show a plots of the thickness of the melted 
layer assuming the rails are made of copper and tantalum, re- 
spectively, for a constant total gun current, which corresponds 
to a unifomly accelerating armature. 

Our model predicts that the thickness of the melted layer 
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Fig. 7. Plot of Eq. (14) , giving the temperature rise, assuming rails made 
from tantalum, plotted vs. position x along the rail. Parameters used are given 
in Tables I and II. The gun velocity is taken as 3000 m/s and the inital 
temperature T = 0°C. Constant current is assumed. 
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Fig. 8. Plot of Eq. JT7}, giving the melt layer thickness vs. x along the rails 
for rails made from copper, tantalum and aluminum. The initial temperature 
of rails is assumed to be T a = 0°C. Constant current is assumed and the gun 
velocity is taken to be 3000 m/s. 



is largest at x = 0, at the breech of the gun rail. This is 
reasonable because the current at the breech of the gun flows 
for the longest time, causing maximum Joule heating at x = 0. 

The intersection of the curves in Figure 6 and 7 with the 
melting point of metal (horizontal line) gives the thickness 
of the melted layer. The surface inside the rails that reaches 
the melting temperature, T me i t , is given by points (x,y) that 
satisfy 

T(x, y, t f ) - T = T^u - T (16) 

Alternatively, we can say that at position x the thickness of the 
rail that reaches the melting point, y me it{x), is given implicitly 



by T(x,y me it,tf) - T 



T„ 



Tq. For any element of 



volume, the temperature rise is due to the length of time that 
the current was flowing through that element. For an element 
at (x, y), current starts to flow at time t(x, y). We assume that 
for all elements the current stops flowing at time tf, when the 
armature leaves the rails. 

From Eq. ( [T6] > we solve for the thickness of the melted 
layer, y m eit(x), as a function of x position along the rail 



Gun Velocity [ J 

sec 

Fig. 9. Plot of Eq.fnJ giving the melted layer thickness for Cu, Ta, and Al 
rails at the gun breech (at x=0) vs. gun velocity, v. The initial temperature of 
rails is assumed to be T a = 0°C. Constant current is assumed. This plot also 
gives an estimate of the melted layer thickness in the armature, see discussion 
in Section VI. 

TABLE I 
Gun parameters. 



quantity 


symbol 


value 




width of rails 


w 


0.10 m 




length of rails (gun length) 


I 


10.0 m 




mass of armature 


m 


20 kg 




derivative of self inductance of rail 


V 


0.46xl0~ 6 


H/m 


magnetic permeability 




4vr X 10- 7 


H/m 



Vmelt{x) = 



8£ 
jJLU V 



1/2 



2 p L' iC (T n 



1 



1/2' 



(I) 

(17) 

Figure [8] shows a plot of the thickness of the melted layer as- 
suming the rails are made of aluminum, copper, and tantalum. 
We used numerical values given in Tables I and II. 

The thickness of the melted layer has a strong dependence 
on gun velocity, which is the velocity of the armature at 



x = £. Figure [9] shows plots of the melted layer thickness 
at x = for Cu, Ta, and Al rails vs. gun velocity. The plots in 
Figure [9] show that for gun velocities above 1500 m/s, using 
materials such a copper, aluminum, and tantalum, the thickness 
of the melted layer in the gun rails increases rapidly with gun 
velocity. 

IV. Current Distribution due to Magnetic Field 
Diffusion 

In this section we obtain the current density in the rail, 
J(x, y, t) (to be used in Eq. (|6]l), by considering the magnetic 
field diffusion into the rails. We treat the diffusion of the 
magnetic field into the rails as diffusion into a plane surface, 
with a time dependent boundary condition on the field, given 
by the armature passing a surface element of the rail. Consider 
an arbitrarily thin armature at position x a at some time t. For 
the geometry in Figure 1, the magnetic field H z is in the 
z -direction. For x > x a , which is outside the rail-armature 
circuit, the magnetic field H z = 0. For < x < x , which 
is inside the rail-armature circuit, the magnetic field has some 
constant value Ho- As the (arbitrarily thin) armature passes the 
point x on the surface of the rail, the surface magnetic field 
changes from H = to some finite value Ho, and the field 
starts to diffuse into the rail. Essentially, as the armature passes 
the point x the boundary condition on the field H z on the rail 
feuWace changes from H z = to H z — Ho- Considering the 

TABLE II 

Material parameters for the rails. 



iclliny tump. ° C conductivity (Ohm m) 



density kg/m 



specific heal J/kg ° C 



[anUilum 



660 
1084 
2996 



3.82 X 10' 
5.8 X 10 7 
7.40X 10 6 



2720. 
8960. 
166(10. 



950. 
440. 
150.62 
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rail to be a plane surface, the magnetic field diffuses into the 
rail according to the equation [12], [17] 



dy 2 



idH z 

hi 



dt 



= 



For a plane surface with magnetic field boundary condition 
H z = for t < and H z = H Q for t > 0, the solution for 

y > is 



iiz = ifn erfc 



2 Jut 



(19) 



where k = l/(p<j) is the magnetic diffusion length, erfc (£) 
is the complimentary error function, erfc (£) = 1 — erf (£) and 
erf (£) is the error function given by 



erf(z) 



dt 



2 

\A ./o 

The current density associated with this magnetic field can be 
obtained from Eq. ([7j, leading to 

,2 



V 'tt nt 



exp 



4«t 



(20) 



where e x is the unit vector in the x-direction, and the auxiliary 
function, J(y,t), is defined by Eq. (20 1. The value of the 
magnetic field is related to the total current I in the rails by 



In 



dy dz J(y, t) = wH 



(21) 



leading to Hq = I /w where w is the width of the rail in the 
z-direction. 

We use the plane surface solution to approximate the 
magnetic field diffusion and resulting current distribution in 
the rail. As the armature sweeps past a surface element in the 
rail at position x, the field starts diffusing into the rail at time 
t (x), where t (x) gives the time the armature passes point 
x. We take the current density to have the form 



J(x, y, z,t) = J (y, t - t a {x)) e x 



(22) 



when the coordinates (x, y, z), and time t satisfy the conditions 

w 

< x < x {t) , y > , \z\ < — , < t < t f (23) 



and otherwise we take 3(x,y,z,t) — 0. Here, x (t) is the 
x-position of the armature at time t, and J(y,t) is given by 



Eq. (20 1. The quantity 3(x, y, z, t)) is an approximation to the 



time-dependent current density in the rails during the EMG 
shot. Note that J(x,y,t) does not depend on z because we 
assume a large (essentially infinite) extent in the z-direction. 

The temperature rise is obtained from Eq. |6]) and we neglect 
the effect of re-distribution of heat during the time of the shot, 
since it takes a long time on the time scale of a shot (which 
is tf), The temperature rise immediately after the shot is 



T ( x iy,tf) -t = 



T O, 



y 



irw 2 pC \ ' 2 tf — t (x) 
where T(a, z) is the incomplete Gamma function, given by 



(24) 



r(a, z) 



t"- 



'd,t 



melting temperature of copper 




X = 1.0 meter 



0.002 0.003 
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Fig. 10. The temperature rise is plotted vs. y coordinate into the rail for 
various positions x along the copper rail. The time dependence of the current 
is determined by the time dependence of magnetic field diffusion into the 
rails. Initial temperature of the rails was assumed to be To =0 C. Parameters 
used are given in Tables [T] and [TT] 



5 0.003 




x position along rail [meter] 



Fig. 1 1 . Plot of the thickness of the melted layer vs. position x along the 
rail, assuming rails made from copper. Curves for three different velocities are 
shown. Initial temperature of the rails was assumed to be To =0 C. Parameters 
are those in Tables III and [TT1 



Equation ( |24[ i is our basic result for the temperature distribu- 
tion in the rails, and is expressed in terms of the current in the 
rails I Q . It is useful to express the temperature rise in terms 
of the gun velocity v. 

In order to get simple results, we assume that the acceler- 
ation is constant, using the relations in Eq. |5]l leads to the 
temperature rise 



T (x>y,tf) - T 



irw 2 pCL'. 



r o, 



fia v 



y 



aVI VI- Jx 



(25) 

Figure 12 shows a plot of the temperature rise given by 



Eq. (25 i, for the case of copper rails, using the parameters 
in Tables I and II. 



In Eq. (25 1, the argument of the r(£) function is large for 

0. For the values 



all values of x and y, except for y near y 
in Tables I and II, the argument of the T function is 



vy 2 fi<j 



54663.71 y 2 
10. - 3.1622VS 



(26) 



Therefore, the r(£) function for £ >> 1 can be approximated 
by its asymptotic expansion 
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(27) 



Near y=0, the behavior of the r(£) function for £ << 1 is 
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Fig. 12. Plot of the thickness of the melted layer vs. position x along the rail, 
assuming rails made from tantalum. Curves for three different velocities are 
shown. Initial temperature of the rails was assumed to be To =0 C. Parameters 
are those in Tables HI and ITU 



logarithmic 



r(o,0 = (-7-£°sE]) + £- 



e 



(28) 



where 7 is Euler's constant, 7 ps 0.5772. 

As before, see Eq. ( |16) , we use Eq.p5|) to solve for the 
thickness of the melted layer in the y-direction by setting the 
temperature rise equal to the melting temperature of the metal. 
Figures 11 and 12 show plots of the thickness of the melted 
layers verses position x along the rail, for several velocities, 
for rails made of copper and tantalum, respectively. From 
these figures, we see that even though Ta has a much higher 
melting temperature than Cu, 2996°C compared to 1084°C, 
respectively, the electrical conductivity of Ta is lower, which 
leads to a comparable thickness for the melted layer for Ta 
and Cu rails. 

V. Armature Heating 

We can estimate the heating of the armature from the 
calculations that we have already done. We now imagine that 
the armature has a finite thickness. Furthermore, depending on 
the EMG design, the armature material can be different than 
that of the rails, and consequently, the electrical conductivity 
a may be different. Unlike the rails, the armature conducts 
current for the whole duration of the shot from t = to t = tf. 
In this sense, the whole length of the armature has current 
flowing in it like the element of rail at coordinate x = 0. The 
thickness of the melted layer of the armature can be found 
from Eq.( 17 1 by setting x = 0, 



ymeit(O) = 



(lav 



1/2 



2 w 2 p l' tc (Tjx 



-To) 



(29) 



Figure 9 shows a plot of the melted layer thickness of the 



armature verses gun velocity, as given by Eq. (29 1 



VI. Efficiency and Melt Layer Thickness 

During an EMG shot, the energy that is supplied to the 
gun appears as projectile kinetic energy, Joule heating of 
the rails and armature, armature and rail heating caused by 
frictional forces between rails and armature, vibration of rails 
and armature, sound and light produced in surrounding air, and 
electromagnetic radiation due to a time-dependent magnetic 



field. We neglect all these effects except for the Joule heating 
of the rails and kinetic energy of the projectile and armature. 



From the current density in Eq. ( 22 1, we calculate the energy 
in one rail, Q, that is dissipated over one shot of the EMG 
by integrating the Joule heating over the time t 1 during which 
the armature is in contact with the rails, 



Q = 2 / d :i x I dt 



■11 
a 



(30) 



where the spatial integration is over all space inside one of the 
rails. The factor of two is due to the fact that there are two 
rails, and by symmetry the energy dissipated is twice that of 
one rail. Assuming constant current, for which the time tf is 
defined in Eq. |5]), assuming the conductivity a is a constant, 
and using the expression given by Eq. (22i for the current 
density derived by considering magnetic field diffusion into 
rails each of length £, we find 
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32 
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(31) 



where we have assumed a constant current (armature accel- 
eration) and used the value of tf in Eq. Q. Note that Q is 
the total amount of energy dissipated in both rails. The two 
forms for Q are obtained assuming the relation between EMG 
velocity and current given in Eq. which assumes a constant 
current. Equation (31 1 gives the amount of energy Q dissipated 



in the rails due to Joule heating over the time tf of one shot. 
Equation (31 1 contains some interesting physics because the 



total energy dissipated as Joule heat is not proportional to 
1^ as is the case for a resistor carrying a constant current 
I Q for a time tf. If the rails acted as a simple resistor with 
resistance R, then the amount of energy dissipated due to 
Joule heating during the shot would be I 2 Rtf. If tf was 
a constant independent of current I a , then the amount of heat 
dissipated in the resistor would be proportional to I 2 . However, 
in an EMG, the time tf is inversely proportional to current 
I , as given by Eq. because the armature spends less 
time in contact with the rails at higher current (due to higher 
average armature velocity at higher current). So if the rails 
acted like a resistor with no skin effect, and the time t f varied 
inversely with current I , then the amount of heat dissipated in 
such a resistor would be proportional to I Q . However, in the 
EMG model that we are considering, we have included the 
skin effect. With increasing current I a , the time tj is shorter 
and there is a stronger skin effect at larger I Q leading to an 
increased resistance R, resulting in an increased Joule heat 
dissipation during the shot time tf, since all current has to 
flow through a thinner channel in the rails. The net effect is 
that the skin effect causes an increased Joule heat over the 
shot, so that for the EMG with skin effect we have Q 



7-3/2 

1 7 



as given in Eq. (31 



The ratio of armature kinetic energy divided by the Joule 
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Fig. 13. The plot shows the thermal q defined in Eq. \32) as a function 
of velocity, using the parameters in Table I, assuming copper rails with 
conductivity in Table II. 



heat is a sort of "thermal q" and is given by 
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The quantity q shows that, in order to achieve a certain velocity 
v, a larger fraction of input energy goes into the armature 
kinetic energy when using shorter rails (smaller £). Figure 15 
shows the dependence of q on velocity v for parameters given 
in Table I and II for gun rails made of copper. However, the 
quantity q does not contain the whole story. We next consider 
the EMG efficiency. 

Within our model, we define the efficiency of an EMG, r\, 
as the armature kinetic energy, divided by the total energy 

1™„,2 



n 



(33) 



\mv 2 + Q + E mag (£) 

We consider the total energy as a sum of three terms: the 
kinetic energy of the armature \mv 2 , the energy dissipated 
by Joule heating Q in Eq. (31) , and the energy stored in the 
magnetic field E mag (£) when the armature is at the end of the 
rails at x = I. In general, the magnetic energy E mag (x) is the 
energy stored in the magnetic field in the rails and surrounding 
space when the armature is at position x: 

1 



1 {x) = -L{x)Ii 



(34) 



As previously defined in Eq. 0, we write L(x) in terms 
of Lo an d L', where Lq is the self inductance of the EMG 
that does not depend on armature position (and corresponds 
to magnetic energy L J 2 /2 stored in the power supply, 
surrounding space, and in electrical leads to the rails) and 
L', which is the derivative of the self inductance with respect 
to armature position x. When the armature reaches the end 
of the rails, the circuit is broken and the magnetic energy 
E mag (£) is dissipated in the form of an electrical arc, sound, 
light, and mechanical vibration. Assuming constant current 
and acceleration, when the armature is at the end of the rails 
at position x = £, the stored magnetic energy is 

\nag{i) = ^ m V 2 (l+ jjj 



E 



(35) 



where we have used Eqs. (4). From Eqs. (30), p3) and ( [35} 
the efficiency of the EMG is then 

1 
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2+Jf 
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(36) 



Equation ( 36 1 shows how the efficiency depends on a number 



of variables. The efficiency of an EMG increases with increas- 
ing electrical conductivity a, velocity v, and rail width w. The 
efficiency depends in a complicated way on the rail length £. 
In the limit of high velocity v, the efficiency of the EMG has 
the limiting value 



lim rj = 



1 



2 4- L " 

u i 



(37) 



The high-velocity limit of EMG efficiency depends on the ratio 
of the stationary self inductance, Lq, to the dynamic part of 
the self inductance of the rails, L' £. We call L' I the dynamic 
part of the self inductance because it depends on the length 
of the rails. For a given gun design, the quantities Lq and L' 
are constants, but in principle the length of the rails, £, can 
be increased to make the term -jfj « 1, thereby increasing 
the high-velocity limit of the efficiency. Longer rails (larger 
£) lead to a higher EMG efficiency at high velocity. However, 
for all velocities the EMG efficiency for this model is always 
less than 1/2. 



The denominator of the efficiency in Eq. ( 36 1 has three 
terms. For the parameters in Table I and II, the third (last) 
term in the denominator is 24.35/-v/(u), where v is in units 
of m/s. Therefore, it is clear that the EMG is more efficient 
at higher velocities. The efficiency will increase with velocity 
significantly when the stationary part of the self inductance, 
Lq, is smaller than the dynamic part of the self inductance, 
L 1 £, so that the term jfj « 1, which may rarely be true 
in real designs unless the length of the rails I can be made 
sufficiently large. 

Stated in another way, a longer EMG is more efficient at 
higher velocities, however, the thickness of the melted layer 
also increases at higher velocities. Therefore, there is a tradeoff 
between efficiency and melting of the rails. Increasing the 
length of the rails may prevent melting, however, this may 
lead to rail length that is not practical for applications. 

In our simple model of an EMG, and our definition of 
efficiency, we have neglected many effects such as friction 
between armature and rails, plasma contacts, and many other 



details. These effects would make the denominator of Eq. (36 1 
larger, and hence would reduce the efficiency. 

VII. Conclusion 

We have constructed a model of the EMG based on the 
electrodynamics of the launch rails and armature. Our model 
takes into account the skin effect in the rails and armature, 
which is significant due to the short-duration and extremely 
high current densities in high-performance EMGs. We used 
two approaches. The first approach, in Section III we modeled 
the skin effect by a constant-in-time current flowing through 
a channel equal to the local skin depth along the rail, 6, 
given by Eq. (7) at each point along the rail. We find that 
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a finite thickness layer of the rails, given by Eq. ( fl7| ), reaches 
the melting point of the metal. Although aluminum is lighter 
by a large factor, we have assumed the rails are made from 
copper because of its significantly higher melting point (1083 
C for copper verses 660 C for aluminum). Figures 6 and 7 
show plots of the temperature rise (immediately after a shot) 
versus depth into the rails, for different positions spaced one 
meter apart along the rails, assuming 10 meter long rails, 
for copper and tantalum rails, respectively. For all positions 
along the rails, immediately after the shot and before thermal 
diffusion takes place, the temperature rise is the highest at 
the rail surface and decreases into the interior of the rails. 
For constant current, near the breech end of the gun the skin 
depth exceeds 15 mm, and a layer approximately 2 mm or 
3 mm thick melts, for tantalum or copper rails, respectively, 
see Figure 8. The model predicts that, immediately after a shot 
and before thermal energy is redistributed, the temperature is 
largest near the rail surface and a thin surface melt layer exists 
for any rail gun, see Eq. ( 17 1. 



In the second approach, in Section IV, we obtained an ap- 
proximate current density by computing the time dependence 
of the diffusion of the magnetic field into the rail surface 
during the EMG shot. From the magnetic field H, we obtain 
the current density from Eq. Q, which is used to compute the 
temperature rise immediately after a EMG shot, see Figures 6, 
7 and 10. Figures 11 and 12 show the melted layer thickness 
we can expect for copper and tantalum rails, for different gun 
velocities (due to different gun currents, which are related to 
velocities by Eq. (4)). We expect this approach (Section IV) 
to be our most accurate estimation of the temperature rise in 
the rails and resultant melting. In both approaches, we find 
the thickness of the melted layer of the rail increases rapidly 
with EMG velocity, see Eq. ( |29] > and Figure 9. 

Finally, in Section V, we estimated the temperature rise of 
the armature, and the resulting thickness of the melted layer. 
The temperature rise in the armature may be a major limiting 
factor in EMG design, see Figure 9. 

In Section VI, we computed the efficiency of an EMG as a 
function of armature velocity and gun length. In this efficiency 
calculation, we only considered Joule heating, kinetic energy 
of the armature, and the stored magnetic energy in the system, 
and we neglected all other energies, such as frictional heating 
between armature and rails. We find that if the stationary part 
of the self inductance is small compared to the dynamic part 
of the self inductance (that of the rails), which is equivalent to 
long rails, the efficiency of an EMG increases with velocity. 
Stated in another way, if the rails are long enough the EMG 
will have an efficiency that increases with velocity. In all cases 
the upper limit of efficiency is 1/2 at high velocity. However, 
with increasing velocity, the thickness of the melted layer in 
the rails also increases, see Figure 9. The thickness of the 
melted layer also decreases with rail length. Therefore, there 
is a tradeoff: for sufficiently long rails, with increasing gun 
velocity the EMG is more efficient, however, at higher velocity 
the thickness of the melted layer is larger, which is likely to 
result in more damage to the rails. Also, we may need rails 
that are too long for practical applications. In our model of 
an EMG, we have neglected many practical effects, such as 



heating and ablation of armature/projectile at high velocities 
in the atmosphere. 

The most important conclusion is that, given the avail- 
able choice of materials for rails and armatures, the gun 
has increased efficiency at higher velocity, however, there is 
increased melting of the rails and damage is more likely during 
each firing of high performance guns, where payload and 
range is designed for (high-velocity) naval guns. This analysis 
was performed for simple rail guns. Other geometries exist, 
although all appear to rely on the sliding armature to carry 
the high currents and hence will generate a characteristic skin 
effect that results in significant rail heating. 

Finally, heating of the rails on a microscopic scale depends 
on inhomogeneities of the metal, such as crystalline alignment 
and dislocations. These inhomogeneities will lead to spatial 
fluctuations of the current density. This will likely lead to 
a spatially inhomogeneous deposit of energy, and local "hot 
spots". Such spatial fluctuations in the case of an EMG are 
likely to lead to the analogue of what is usually termed "gun 
barrel erosion". A comment made at the end of the chapter 
written by I. Ahmad in 1988, entitled "The Problem of Gun 
Barrel Erosion: An Overview", is still highly relevant: "Other 
forms of erosion/corrosion problems might appear as a result 
of advances made in the development of liquid propellant guns 
and electromagnetic gun technologies. In fact, the erosion in 
the latter could be quite severe, as it involves interaction of 
high-temperature plasma with the launch surface causing it 
to partially melt at each firing event. It will be necessary to 
identify materials and design methodologies to minimize these 
problems." Our paper does not include any effect of the high- 
temperature plasma, but it shows the effect of Joule heating 
and geometric design considerations on the rails. 
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